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15 Abstract: 

16 In the present work, we determine how three well documented models of the 

17 magnetospheric electric field, and two different mechanisms proposed for the formation 

18 of the plasmapause influence the radial distance, the shape and the evolution of the 

19 plasmapause during the geomagnetic storms of 28 October 2001 and of 17 April 2002. 

20 The convection electric field models considered are: Mcllwain’s E5D electric field 

21 model, Volland-Stem’s model and Weimer’s statistical model compiled from low-Earth 

22 orbit satellite data. The mechanisms for the formation of the plasmapause to be tested are: 

23 (i) the MHD theory where the plasmapause should correspond to the last-closed-. 

24 equipotential (LCE) or last-closed-streamline (LCS), if the E-field distribution is 

25 stationary or time-dependent respectively; (ii) the interchange mechanism where the 

26 plasmapause corresponds to streamlines tangent to a Zero-Parrallel-Force surface where 

27 the field-aligned plasma distribution becomes convectively unstable during enhancements 

28 of the E-field intensity in the nightside local time sector. 

29 The results of the different time dependent simulations are compared with concomitant 

30 EUV observations when available. The plasmatails or plumes observed after both 

31 selected geomagnetic storms are predicted in all simulations and for all E-field models. 

32 However, their shapes are quite different depending on the E-field models and the 

33 mechanisms that are used. Despite the partial success of the simulations to reproduce 

34 plumes during magnetic storms and substorms, there remains a long way to go before the 

35 detailed structures observed in the EUV observations during periods of geomagnetic 

36 activity can be accounted for very precisely by the existing E-field models. Furthermore, 

37 it cannot be excluded that the mechanisms currently identified to explain the formation of 

38 “Carpenter’s knee” during substorm events, will have to be revised or complemented in 

39 the cases of geomagnetic storms. 

40 

41 

42 1. Introduction. 

43 

44 The plasmapause position prediction is very important for the scientific community 

45 because its location determines the influence of wave-particle interaction processes on 
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1 radiation belts (RB) formation. Cornwall et al. [1971] noted that maximum instability of 

2 ElectroMagnetic Ion Cyclotron (EMIC) waves generated by ring current (RC) anisotropy 

3 should occur just within the plasmapause. He also suggested that pitch-angle diffusion of 

4 RC ions resonant with the EMIC waves generated within the plasmapause should be " ,J 

5 important loss process for the RC. Recent global self-consistent RC/EMIC waves 

6 modeling by Khazanov et al. [2006] conforms this initial study. It was found that global 

7 EMIC wave distribution has a highly plasmapause organized structure during the May, 

8 1998 storm period that has been presented in this study. Khazanov et al. [2007] also 

9 found that plasmaspheric energy deposition to the thermal electrons via Landau EMIC 

10 wave damping is very critical to position of plasmapause. 

11 

12 During magnetic storms, under certain conditions, relativistic electrons with energies >1 

13 MeV can be removed from the outer RB by EMIC wave scattering. Recent calculations 

14 suggest that pitch angle scattering via EMIC waves can compete with Dst effect as a 

1 5 mechanism for depleting relativistic electrons from the outer RB zone during the initial 

16 and main phases of a magnetic storm [Summers and Thome, 2003; Albert, 2003]. So, it 

17 becomes more and more obvious that RC and RB populations are very sensitive to the 

1 8 core plasmasphere distribution and specifically to the position of the plasmapause 

19 [Spasojevic et al., 2004; Baker et al., 2004], 

20 

2 1 The position of the plasmapause is determined by the interplay between the co-rotation 

22 and convection electric fields. This statement is very general but it does not really 

23 identify and grasp the specific physical mechanism that is underlying the formation of the 

24 plasmapause. The magnetospheric convection electric field, controlled by the solar wind 

25 conditions and the level of geomagnetic activity, is a key factor in all existing theories for 

26 the formation of the plasmapause. Therefore, it is important to have first a reliable 

27 magnetospheric electric field model. Since there is no way to determine directly, at each 

28 instant of time, the global electric field distribution in the whole magnetosphere, various 

29 empirical and mathematical models have been built, tentatively, with various degrees of 

30 sophistication. The E-field models used in the present paper are (i) the VSMC model 

3 1 introduced by Volland and Stem [Volland, 1973; Stem, 1975] and adapted by Maynard 
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1 and Chen [1975], (ii) the E5D model derived by Mcllwain [1986] from ATS5 and ATS6 

2 observations at geosynchronous orbit, and (iii) the model of Weimer [1996] determined 

3 from ionospheric measurements; these empirical electric field models are respectively 

4 associated with (i) a centered dipole magnetic field, (ii) the M2 magnetic field determined 

5 from geosynchroneous measurements, and (iii) Tsyganenko’s [1996] magnetic field 

6 model determined from various statistical magnetometric in-situ measurements. These 

7 models will be briefly described in the next section. 

8 In addition to a reliable magnetospheric electric field model, one needs also to have a 

9 correct physical theory for the formation of the plasmapause, i.e. a specific physical 

10 mechanism that accounts for the observations. Several such mechanisms have been 

1 1 proposed in the past and can be simulated numerically once an E-field model has been 

12 adopted. 

13 

14 In the following sections we use the three E-field models mentioned above, and compare 

1 5 their influence on the dynamics of the plasmasphere during the geomagnetic storm on 28 

16 October 2001. Moreover, for the E5D model, the dynamical evolution of the cold plasma 

17 is tracked by using two different categories of simulations to predict the positions of the 

1 8 plasmapause and its deformations at different times precisely when EUV observations are 

19 available. 

20 

21 There are other sophisticated models for the magnetospheric electric field distribution: 

22 for instance, AM3E developed by Richmond and Kamide [1988]. This popular E-field 

23 model was obtained from ground-based (magnetometer and radar) and ionospheric 

24 (DMSP satellite) data. Although it might be of some interest to use this additional model, 

25 it is not completely clear that ionospheric electric field distribution can be directly 

26 mapped up into the equatorial region of the magnetosphere where plasmapause knee start 

27 to form. Not only is the actual 3D distribution of magnetic field lines not completely 

28 guaranteed, but the existence of field-aligned electric potential drops is likely to 

29 jeopardize such a field-aligned mapping of ionospheric E-field into the magnetosphere. 

30 
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1 In several studies, different convection electric field models have already been compared 

2 with observations [Jordanova et al., 2001; Boonsiriseth et al., 2001, Khazanov et al., 

3 2004]. Their influence on the ring current was analyzed in detail by Chen et al. [2003]. 

4 But their effects on particles of kftver energies (<1 0 eV) populating the plasmasphere was 

5 less well studied, except by Liemohn et al. [2004] who compared some electric field 

6 models, and their effects on plasmaspheric morphologies; he used (i) a modified 

7 Mcllwain electric field model, (ii) the Weimer model, (iii) as well as another self- 

8 consistent electric potential model for the time span of the recovery phase of the 17 April 

9 2002 magnetic storm. These authors found that all these models have certain strengths 

10 but also some weaknesses in predicting the observed plasmapause position during this 

1 1 storm. They found especially that the electric field intensity of Weimer’ s model was a bit 

12 too strong in the inner magnetosphere, leading to a too small plasmasphere. Liemohn’s 

13 modified Mcllwain model (which differs from the original E5D model, and was not used 

14 with the associated M2 magnetic field model) has a too small electric field intensity 

15 around noon, leading to a plasmapause position that does not correspond to the EUV 

16 observations on the dayside although a good fit was obtained in the nightside. 

17 

18 In the present work, we also present the results obtained during this magnetic storm of 1 7 

19 April 2002, but we have focused our attention on another clean and isolated geomagnetic 

20 storm: that observed on 28 October 2001 . The VSMC, unmodified E5D and Weimer E- 

21 field models have been used with the hope to determine which of them might be the most 

22 appropriate for this magnetic storm. 

23 

24 It is commonly assumed that the topologies of both the electric field and magnetic field 

25 determine the position where the plasmapause is formed at the time of geomagnetic 

26 storms and substorms. Postulating that the plasmapause coincides with the last closed 

27 equipotential of the magnetospheric electrostatic field distribution as proposed by Brice 

28 [1967], this led authors to derive an electric field topology from observed plasmapause 

29 positions: e.g. Maynard and Chen [1975]. This led also Goldstein et al. [2002] to add an 

30 electric field component penetrating into the plasmasphere to obtain plasmaspheric 

31 shapes and positions fitting those observed on 24 May 2000, when a plasmaspheric 
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1 shoulder was found in the EUV observations. In a subsequent study, Goldstein et al. 

2 [2003] included an additional E-field distribution related to sub-auroral polarization 

3 stream to obtain a better fit with plasmapause positions observed by EUV/IMAGE on 2 

4 June 2001. 

5 The MHD simulations are essentially based on the postulate that the plasmapause 

6 coincides with the Last Closed Equipotential (LCE) of a stationary global 

7 magnetospheric electric field. When the magnetospheric electric field is not stationary, 

8 the plasmapause can not be identified with the LCE of the time-dependent E-field, but 

9 may be assumed to correspond to a Last Closed Streamline (LCS). This is what will be 

10 assumed in the MHD simulations presented below. The second mechanism used below to 

1 1 predict the plasmapause positions and shapes is based on the mechanism of interchange 

12 reported in the book by Lemaire and Gringauz [1998], These mechanisms are recalled 

13 and described in section 3. 

14 

15 In section 4, we illustrate the plasmapause positions simulated for all three different 

16 electric fields models and with both kinds of mechanisms. The results are then compared 

1 7 with observations of EUV/IMAGE. Discussion and conclusion are given in the last 

1 8 section. 

19 

20 2, Models of electric and magnetic fields 

21 

22 2.1 Volland-Stern’s and Maynard-Chen’s convection electric field (VSMC) 

23 

24 The Volland-Stem model [Volland, 1973, Stem, 1975] is a simple mathematical model 

25 where a uniform dawn-dusk convection electric potential distribution is applied across 

26 the magnetosphere. It has become very popular because of its simplicity and portability. 

27 In this model there is no induced electric field resulting from time dependent magnetic 

28 field variations as can be envisaged during strong geomagnetic storms; the 

29 magnetospheric electric field derives from a scalar potential which, in a co-rotating frame 

30 of reference, is given by 

31 O = AR 2 sin^ 
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where R is the equatorial radial distance, (|) is the azimuthal angle from noon and the Kp 


dependent factor A = 


0.045 


(l-0.159K/? + 0.0093Kp 2 ) 3 


kV_ 

Rl 


determines the convection 


electric field intensity. 


The Kp dependence of this empirical model was obtained by Maynard and Chen [1975] 
by adjusting the last closed equipotential (LCE) of the total E-field with plasmapause 
positions determined by OG03 and OG05 satellite observations; this is why this model 
has been given the acronym VSMC in the following. The Kp dependence is very 
sensitive especially at small values of Kp; the value of A varies from A=45 V/ R E 2 for Kp 
= 0 to over 800 V/Re 2 for Kp = 6. 


Fig. 1 illustrates the equatorial contour maps of the equipotentials for the VSMC 
convection electric field every two hours from 0:00 UT to 10:00 UT during the 
geomagnetic storm of 28-10-2001. It shows how the dawn-dusk electric field component 
is shielded close to the Earth and how it is assumed to change over this lOh period. 
During this period of time, Kp increased from V at 0:00 UT to T at 05:00 UT, as 
illustrated in the top panels of Fig. 9. 

In our simulations, this convection electric field is used with a centered dipolar magnetic 


field: B = 


3100 

R 2 


. Indeed, this was the magnetic field assumed by Maynard and Chen 


— > ■ — > 

[1975] to compute the MHD convection velocity: Ex B! B 2 . 


2. 2. Mcllwain’s E5D convection electric field 


Another analytical representation of the magnetospheric convection electric potential was 
derived by Mcllwain [1986] from electron and proton dynamical spectra measured at 
geosynchronous orbit during the ATS-5 and ATS-6 missions: 
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2 with R ar = 9. 8 - 1. 4 cos 0 + (-0.9 -0.3 cos ^)- — — — 

ar r l + OAKp 

3 

4 The E5D model depends also on the three-hourly geomagnetic activity index Kp. It has 

5 been found in past statistical studies that this index controls the observed positions of the 

6 plasmapause, when the level of geomagnetic activity changes. Note however that, unlike 

7 the VSMC model, the E5D model was not derived by adjusting its parameters to fit 

8 observed plasmapause positions, but was directly deduced from ATS-5 and ATS-6 

9 particle flux measurements at geosynchronous altitude, precisely in the region where the 

10 plasmapause is formed. The constants in the E5D electric potential model have been 

1 1 adjusted to fit the dynamical energy spectra in the range from 1 keV to 100 keV of 

12 electrons and protons injected following substorm events and observed in the equatorial 

13 region; according to Mcllwain, this model is not necessarily reliable for Kp > 6, nor in 

14 the case of geomagnetic storms when rapidly changing electric field intensities are 

1 5 induced inside the magnetosphere. 

16 It should be emphasized that Mcllwain [1986] did not develop his E5D model based on 

17 ATS-5 and ATS-6 observations collected during geomagnetic storm events, when Dst has 

18 large excursions. In principle, this model is designed to simulate the plasmapause 

19 formation and deformations during relatively steady state conditions following substorm 

20 injection events. Nevertheless we will use it here during a geomagnetic storm with the 

21 caution that if the plasmapause positions predicted with this model would not fit 

22 adequately the observations during the geomagnetic storm, this may be due to the 

23 inadequacy of the E5D model to fit the actual electric field distribution during this period 

24 of time, and not necessarily to the inadequacy of the physical mechanism assumed to 

25 form the plasmapause. 

26 
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1 Fig. 2 illustrates the equatorial equipotential contours for the E5D convection electric 

2 field model. These maps are shown every two hours of UT from 0:00 to 10:00 during the 

3 geomagnetic storm of 28-10-2001. It clearly indicates the shielding of the dawn-dusk 

4 electric field component near the Earth in the dawn sector. A comparison with Fig. 1 

5 indicates that the E5D electric field model is less sensitive to changes of Kp than the 

6 VSMC model. 

7 To force a stronger Kp-variation, Liemohn et al [2001] re-scaled the E5D model to match 

8 the cross-polar cap potential difference to other types of observations. This led them to 

9 define a “modified Mcllwain” E-field model that should no more support ATS-5 and 

10 ATS-6 observations. In our simulations, we keep the original E5D model, as well as the 

1 1 associated M2 equatorial magnetic field model since both were derived from concomitant 

12 ATS-5 and ATS-6 observations. 

13 This M2 magnetic field model has a day-night asymmetry, and must be used in 

14 association with E5D to account for the actual day-night asymmetry of the convection 

15 velocities. Figure 3 illustrates the iso-contours of the equatorial magnetic field intensity 

16 corresponding to the M2 model. Beyond 5 earth radii, the deviations from a s imp le 

17 symmetric dipole become large especially in the nightside sector of the magnetosphere. 

18 Note that this B-field model does not depend on Dst nor on Kp unlike the model of 

19 Tsyganenko [1996], Of course, this may be viewed as another limiting factor 

20 compromising its application to study the magnetic storm of October 28, 2001 . 

21 

22 2.3. Weimer’s convection electric field 

23 

24 Unlike the two previous E-field models, Weimer’s model does not depend on the 

25 geomagnetic activity level Kp. Weimer's [1996] E-field model is driven by solar wind 

26 parameters: interplanetary magnetic field magnitude, solar wind velocity and dipole tilt 

27 angle. It was derived from low altitude ionospheric convection velocity measurements at 

28 high latitudes. Unfortunately, these observations are collected far from the equatorial 

29 region where the plasmapause is formed during substorm events. 
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1 Weimer electric potential [Weimer, 1996] is given by an expansion in spherical 

3 ) 

2 harmonics: <&(#, <p) = (A lm cos m<p + B lm sin m (p)P? (cos 9) 

1=0 

3 where 0 is a function of the geomagnetic co-latitude,' (j) is the magnetic local time (MLT) 

4 and P™ are the associated Legendre functions. The coefficients Ai m and Bi m were derived 

5 by a least error fit lfom multiple satellite measurements of the ionospheric convection 

6 velocity. 

7 

8 Fig. 4 illustrates the equatorial contour maps of Weimer’s convection electric potential 

9 every UT hour from 00:00 UT up to 1 1 :00 UT during the geomagnetic event of 28-10- 

10 2001. These equatorial isocontours of Weimer’s equipotential are quite different from 

1 1 those illustrated in Figs. 1 and 2. The shielding is less efficient in the dawn sector than at 

12 dusk unlike in Fig. 2. From the distance between equipotential lines, it can also be seen 

13 that the electric field intensity is stronger than those illustrated in Figs. 1 and 2. 

14 Moreover, Weimer’s model depends on solar wind parameters varying over smaller time 

1 5 scales (60 minutes) and with larger amplitudes than the three-hourly index Kp controlling 

16 the VSMC and E5D models. Weimer’s electric field intensity becomes especially strong 

17 in the dusk MLT sector at 4:00 UT at the beginning of the geomagnetic event while for 

1 8 the E5D model it is strongest in post-midnight sector (Fig. 2) and in both the dawn and 

19 dusk sectors for the VSMC model (Fig. 1). 

20 

21 Weimer’s electric field is used in association with Tsyganenko’s [1996] magnetic field 

22 model that is controlled by the following input parameters: solar wind pressure, Dst, the 

23 Y and Z components of the interplanetary magnetic field, and the geodipole tilt angle. 

24 

25 2.4. Co-rotation electric field 

26 

27 To determine the total magnetospheric electric field in a non co-rotation frame of 

28 reference, the co-rotation electric field must be added to the convection electric field, 

29 postulating that these two E-field distributions could be superposed as in vacuum, despite 

30 the fact that the dielectric constant and permittivity of plasmas is much larger than that of 
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1 free space. However, according to Vasyliunas [2001] this assumption is questionable, 

2 since any (external) electric field imposed from outside of a plasma system does directly 

3 determine the convection velocity deep inside this plasma system. It is the plasma bulk 

4 motion and the generalized Ohm’s law that determine the E-field inside the 9 

5 magnetosphere; this internal E-field distribution does not necessarily coincide with a 

6 ‘simple’ superposition of the ionospheric co-rotation electric field, and an external 

7 convection electric field induced by the solar wind. Indeed, due to the large dielectric 

8 constant, the latter does not penetrate inside the magnetosphere. Although such a simple 

9 superposition in a plasma (i.e. a highly dielectric medium) is not consistent with classical 

10 electrodynamics of dielectric material, it seems, nevertheless, to approximate the actual 

1 1 magnetospheric E-field distribution with mitigated success, at least this is what was 

12 generally considered within the community. 

13 

1 4 The potential of the co-rotation electric field is given by: 

15 <& = -—[kv\ . 

R 

16 

17 Note that it has been inferred from IMAGE observations that the azimuthal velocity of 

1 8 the plasmasphere is often slower (10% in average) than co-rotation in the outermost 

19 layers of the plasmasphere [Burch et al., 2004], This effect has been tentatively taken into 

20 account in some simulations by reducing the co-rotation potential by this ad hoc factor 

21 [Pierrard, 2006], In the present work, however, we will consider that co-rotation is 

22 applicable, since all three E-field models have been derived, originally, under such an 

23 assumption. 

24 In any case, the co-rotation electric field dominates near Earth, where the equipotential 

25 lines are closed and almost circular. This is evidenced in Figs. 5, 6 and 7 showing the 

26 equatorial contour maps of the total electric potential in a non-corotating frame of 

27 reference, respectively for VSMC, E5D and Weimer’s convection electric field models 

28 every two hours on 28 October 2001 . 

29 
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1 Note the significant differences between all three models. The last closed equipotential 

2 has a stagnation point at 1 8:00 MLT in the dusk sector for VSMC and E5D models, while 

3 it is located at later MLT, in the post-dusk local time sector, for Weimer’s model. The 

4 LCE is everywhere closer to the Earth for the VSMC model than for the E5D one, for any 

5 value of the geomagnetic activity level Kp > 1 . 

6 The E-field intensity increases to larger value for the Weimer model than for the two 

7 other models during this geomagnetic storm. Therefore, the last closed equipotential 

8 penetrates closer to Earth for Weimer’s model than for the two other models. For the 

9 VSMC model, the maximum E-field intensity is located in the morning sector (at 06:00 

10 MLT). In the E5D model, the maximum intensity is in the post-midnight sector (around 

1 1 2:00 MLT). Weimer’s model also shows a maximum intensity in the post-midnight sector 

12 during the first hours preceding the geomagnetic storm similar to that displayed in the 

13 E5D model. During the storm main phase, this peak value moves closer to 06:00 MLT in 

14 both models. 

15 

16 The evolution of the electric field (and magnetic field) in Weimer’s model is controlled 

17 by the solar wind and geomagnetic parameters corresponding to 28 October 2001. Fig. 8 

18 shows the variation of the solar wind density of protons, the solar wind velocity and the 

19 Bx and By components of the interplanetary magnetic field obtained from SPENVIS 

20 (www.spenvis.oma.be) from 27 to 30 October 2001 . The three hourly Kp, hourly Dst 

2 1 and Bz values are given in the three top panels of Figs. 9 to 14. 

22 Note that the 28 October 2001 storm is clearly associated to an increase of the solar wind 

23 velocity and density. Moreover, like most geomagnetic storms it is initiated by a 

24 southward turning of the interplanetary magnetic field direction. The geomagnetic 

25 activity level Kp index increases up to T. The Dst index decreases by more than 150 nT. 

26 

27 3. Dynamical simulations 

28 

29 3. 1. Ideal MHD simulations 

30 
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1 Different mechanisms have been proposed to form plasmaspheric “knee” or plasmapause. 

2 These mechanisms are described in detail in the book by Lemaire and Gringauz [1998]. 

3 The early theoretical MHD simulations by Grebowsky [1970] predicted that the evolution 

4 of the plasmapause is determined by the? ideal MHD motion of the LCE at an arbitrarily 

5 chosen initial time to. In Grebowsky’ s early dynamical simulations, the plasmapause was 

6 assumed to coincide with the LCE surface at to. But at any subsequent instant of time, the 

7 plasmapause did not coincide with the LCE of the changing convection E-field, unless 

8 the magnetospheric convection electric field would be independent of time. 

9 Our ideal MHD simulations differ from those of Grebowsky. They resemble more closely to 

10 those developed by Rasmussen [1992]. We launch plasma elements at 23:00 MLT, every 

11 0. 1 5 Re at radial distances ranging from 1 .2 to 6 Re along an equatorial radius. The drift 

-4 — > t ) 

12 path of these ideal MHD plasma elements are calculated from Ex B! B 2 where E and B are 

1 3 respectively the electric field and the magnetic field given by the adopted models. We let the 

14 plasma drift around the Earth for 24 h, and determine the last closed streamline (LCS) which 

15 is then used as the initial plasmapause position at time t G . The ideal MHD drift path of all 

1 6 plasma elements proceeds until a time when EUV observations are available and can be used 

17 for comparison. Plasma elements are continuously launched from 23:00 MLT every 5 

1 8 minutes, and tracked until they are deviated near the stagnation point, and are eventually lost 

19 at the magnetopause; we stop tracking any plasma element once it has completed one turn 

20 around the Earth. The equation of motion of each plasma element is integrated by a Runge- 

21 Kutta method with an adapted time step to satisfy a correct numerical accuracy; the 

22 equatorial positions of all plasma elements are stored every 5 min to map the deformations of 

23 the LCS as a result of Kp variations. For these ideal MHD simulations the LCS is assumed to 

24 coincide with the plasmapause; it should be reminded that the LCS does not coincide at any 

25 time with the LCE of the convection electric field. 

26 

27 In the regions close to the Earth where co-rotation is enforced, the plasmaspheric plasma is 

28 always highly coupled with that of the ionosphere. Since all flux tubes are on closed 

29 streamlines in this region, they are assumed to be completely full, and their densities are 

30 given by the empirical model of Carpenter and Anderson [1992]: 


12 



1 log 10 N = -0.3 145L + 3.9043 where N eq is the number density in electrons/cm 3 in the 

2 equatorial plane; it is independent of the local time angle, but it is a function of the Mcllwain 

3 parameter L. Note also that it does not depend on the level of geomagnetic activity. 

4 The plasma elements located further away from the Earth are more significantly affected by 

5 the solar wind induced magnetospheric convection. When the convection electric field is 

6 enhanced during periods of substorm injection events, the core of the plasmasphere where 

7 co-rotation dominates shrinks. The outer flux tubes beyond the last closed streamline (LCS) 

8 are eventually depleted in less than 24 hours. When the convection electric field intensity 

9 decreases, the outer flux tubes can gradually refill with cold plasma flowing up from the 

10 ionosphere. 

11 

12 Figs. 9 (three upper panels) illustrate the results of these MHD simulations for 28 October 

13 2001 at 00:00 UT for the three electric field models described above. In the all following 

14 Figures, the upper left panels show the results for the E5D model; upper middle panels those 

15 for VSMC model, and upper right panels for Weimer’s model. The symbols represent the 

1 6 plasma elements corresponding to the flux tube contents that were launched at 23 :00 MLT 24 

17 hours before to. At small radial distances, the symbols form concentric circles corresponding 

18 to a contour of constant density content. At large radial distances, the contours of constant 

19 tube content become elongated in the pre-dusk sector, indicating a decrease of the number 

20 density in this sector. Bottom left panel in Fig. 9 corresponds to the result of the simulation 

21 based on the interchange mechanism for the formation of the plasmapause and E5D E-field 

22 model. 

23 

24 We compare the successive positions of the plasmapause found with the MHD 

— — ► 

25 simulations assuming each point moves with the instantaneous Ex B/ B 2 drift velocity 

26 calculated with the three different E-field models: VSMC, E5D, and Weimer. Note that 

27 this convection velocity is never parallel to the equipotential lines of these models unless 

28 the E-field and B-field would both be stationary. Due to the continuously changing E- 

29 field distribution, the position of the LCS is a function of time. Let us recall that the LCS 

30 corresponds to the envelope of all drift paths that have not reached the magnetopause 
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1 boundary during the previous 24 horn's. This is why it should be denoted LCS-24h. A 

2 different LCS could have been obtained if the closure time would have been 25 hours, 2 

3 days, 3 or 6 days as in the MHD simulations of Chen and Wolf [ 1 972]. 

4 The position of the LCS-24h found" with our MHD simulations corresponds to the boundary 

5 where the number of plasma elements per unit area decreases sharply in Figs 9 to 15 (except 

6 at 23:00 MLT and beyond where new plasma elements are continuously launched every five 

7 minutes). If the plasmapause is identified with one of the last closed streamlines, it should be 

8 specified what is the closure time that has been adopted in the MHD simulations. This means 

9 that the location of the plasmapause, a physical boundary, would depend on the arbitrary 

10 choice of a closure time in any MHD simulation. This is clearly a conceptual limitation of the 

1 1 plasmapause with the LCS and even worse with the LCE. 

12 

13 Figs. 9 to 14 show that a plume is produced in the dusk sector during the magnetic storm in 

14 all three MHD simulations. This plume in the LCS-24h is similar to that observed by 

15 EUV/IMAGE after 18:00 UT. 

16 Fig. 10 shows that at 04:30 UT, the LCS-24h has been pushed inwards in the nightside 

17 region. Although the EUV observations show also an inward shift of the plasmapause near 

1 8 midnight, the overall shape of the plasmasphere is however very different from that of the 

19 LCS-24h obtained with the MHD simulations. While the LCS has rather smooth and regular 

20 shape, the actual plasmapause is much more indented and irregular. This can be a 

21 consequence of the rather poor time resolution of the Kp index controlling the E-field 

22 distributions. Since the magnetospheric E-field is changing over smaller time scale, the actual 

23 MLT variation of the LCS is expected to be more irregular than illustrated in Figs. 10 to 14. 

24 Another important consequence of the shorter time variation of the E-field distribution, and 

25 of the associated plasma flow which must be quite variable, was pointed out by Dungey 

26 [1967] soon after the discovery of the sharpness of “Carpenter’s knee”. He wrote : “ Some 

27 tubes of high density (i.e. from inside the plasmasphere) should sometimes be swept out on 

28 the day side, and some tubes, after entering from the tail, should enter the inner region and, 

29 after a few days, should have intermediate values of density. It then seems rather surprising 

30 that the knee should be so sharp, but the variable model would predict a patchy density in the 

3 1 region near the knee and this could be the true state” (sic). We share Dungey’s early concern 
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1 about the sharpness of a knee along the LCS. The short time scale variability of the 

2 magnetospheric plasma flow does not support the formation of sharp density gradients like 

3 those observed soon after the peeling off of the plasmasphere which are associated with 

4 substorm events or inward motion of iiij ection boundaries. However, the variability of the 

5 magnetospheric plasma flow may well account for the existence of the Plasmasphere 

6 Boundary Layer where Carpenter and Lemaire [2004] identified patchy plasma density 

7 irregularities. 

8 

9 Note that the LCS corresponds in the MHD simulations of Chen and Wolf [1972] to the 

1 0 limit between magnetic flux tubes which have refilled for more than 24 hours and those 

1 1 that have refilled for less than 24 hours. If such a limit should correspond to the actual 

12 plasmapause, it must be admitted that the position of this physical boundary would 

13 depend on the choice of the refilling time which in our simulations was arbitrarily taken 

14 to be 24 hours. In Chen and Wolfs MHD simulations different plasmapause locations 

15 were obtained when the refilling time was assumed to be 6 days, 5 days, 4 days, 3 days, 2 

16 days or 24 hours. 

17 This illustrates again the drastic limitation of identifying the plasmapause by any kind of 

1 8 MHD simulations whose results depend on an arbitrary choice of the closure time and a 

19 refilling time. Furthermore, the results of these MHD simulations rely also on the 

20 arbitrary choice of the initial time, to, when such simulations are assumed to start, and on 

21 the assumption of the density distributions in all flux tubes at this starting time. 

22 

23 3.2. Interchange mechanism 

24 Another mechanism has also been proposed for the formation of the plasmapause: the 

25 interchange mechanism, illustrated by Lemaire and Kowalkowski [1981]. According to 

26 this mechanism, the plasma becomes unstable above a Zero Parallel Force (ZPF) surface 

27 where the parallel component of the gravitational plus centrifugal accelerations is equal 

28 to zero [Lemaire, 1985, 2001]. This occurs in the nightside local time sector during 

29 substorms and storms at an equatorial distance that depends on the distribution of the 

30 convection velocity. The plasmasphere is peeled off in this sector by the centrifugal effect 

3 1 which is enhanced when the convection velocity is enhanced beyond the inner edge of a 
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1 substorm injection boundary. By this mechanism, a new plasmapause forms closer to the 

2 Earth than the LCS determined by the MHD simulations presented and discussed above. 

3 This mechanism has been used in numerical simulations performed by Lemaire [1985, 

4 2000] and by Pierrard and Lemaire [2004], The positions of the plasmapause predicted by 

5 this mechanism with the E5D model have been compared to observations of IMAGE for 

6 typical dates [Pierrard and Cabrera, 2005, 2006; Pierrard, 2006]. 

7 

8 In these non-MHD simulations based on the interchange mechanism, plasma holes with a 

9 density smaller than the background density are launched at 23:00 MLT in the equatorial 

1 0 plane. Due to the interchange motion, these holes drift ultimately toward an asymptotic 

1 1 trajectory where the radial component of the gravitational force and centrifugal force 

12 balance each other. This corresponds to the Zero Radial Force (ZRF) surface. However, 

13 the parallel components of these forces balance each other closer to the Earth along a 

14 virtual surface that Lemaire [1985] has called the Zero Parallel Force (ZPF) surface. It is 

1 5 assumed that the field aligned plasma distribution becomes convectively unstable along 

1 6 all flux tubes that traverse or are tangent to this virtual ZPF surface. In these flux tubes 

17 field aligned flow velocity is enhanced, and the plasma density is reduced due to its 

1 8 upward expansion. For a dipole magnetic field distribution the equatorial distance of the 

19 ZPF is 3 2 3 time smaller than that of the ZRF surface. We assume that for the M2 

20 magnetic field model the minimum equatorial distance of the ZPF surface is also 

21 approximately 3 2/3 times smaller than that of the ZRF surface, as it is for a dipole B -field. 

22 

23 At time to we start with a plasmapause whose position has been determined at all MLT 

24 angles by the history of the geomagnetic activity level during the previous 24 h. Note that 

25 the simulations could also be started by using an observed plasmapause position at the 

26 time to like in the study of Goldstein et al. [2003]. This procedure would of course favor 

27 the agreement between the results of the simulation and the observations at any 

28 subsequent times. Unfortunately, this procedure does not tell us about the physical 

29 mechanism that has formed this initial plasmapause; it may only be useful to test the 

30 appropriateness of the E-field models chosen in simulations, after the plasmapause had 

3 1 already formed. 
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1 

2 In the next section, we describe and compare the results obtained with the following 

3 series of simulations: MHD with the E5D model (MHD E5D), MHD with the Volland- 

4’ Stem-Maynard-Chen model (MHD VSMC), MHD with the Weimer model (MHD W), 

5 and interchange with E5D model (IC E5D). 

6 

7 4. Discussion of the results 

8 

9 4.1 Magnetic storm of 28-10-2001 

10 Figs. 9 to 14 illustrate the equatorial position of the plasmapause obtained with two 

1 1 different mechanisms for the formation of the plasmapause and three different electric 

12 field models, respectively at 0:00, 4:30, 8:00, 12:00, 1 8:00 and 20:00 UT. The 

13 simulations correspond to MHD E5D (left upper panels), MHD VSMC (middle upper 

14 panels) and MHD W (right upper panels). The results obtained with the interchange 

15 mechanism and the E5D model are also presented in the left bottom panel d of these 

16 graphs. The plasmapause found with the interchange mechanism is closer to Earth than 

17 with MHD simulations when a same E-field model is used. We do not provide 

1 8 simulations with interchange mechanism for the other electric field models since this 

19 would prohibitively increase the size of this paper without improving significantly its 

20 main results. Note also that the use of the interchange mechanism would not be 

21 physically significant in association with an E-field like the VSMC model which has 

22 been fitted in order to match the LCE with observed plasmapause positions. 

23 

24 The top panels of Fig. 9 to 14 illustrate the Kp and Dst indices and the northward 

25 component of the interplanetary magnetic field (Bz) from the beginning of 27-10-2001 to 

26 the end of 29-10-2001. At 00:00 UT, the level of geomagnetic activity is very low (Kp = 

27 1) and increases to 7" over the next 3 hours when the main phase of the geomagnetic 

28 storm starts, i.e., when Dst begins a gradual drop, as shown on Fig. 9. 

29 

30 At 04:39, 08:03, 18:05 and 19:58 UT, there are exploitable observations from the EUV 

3 1 in s trument since the orbit of the satellite IMAGE is close to its apogee. These images are. 
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1 presented in the bottom right panels of Figs. 10, 11, 13 and 14, respectively. These 

2 observations are intensity maps of the 30,4 nm emissions of Helium ions integrated along 

3 the line of sight. They are projected in the geomagnetic equatorial plane in the SM 

4 coordinate system with the program XForm available at 

5 (ftp://euv.lpl.arizona.edu/pub/bavaro/unsupported/') . This software tool enables to view 

6 the plasmapause cross section from over the North Pole like in our simulations. The 

7 plasmapause is assumed to be the sharp edge where the brightness of 30.4 nm He + 

8 emissions drops drastically. This boundary is illustrated by a white line corresponding to 

9 a threshold equal to 40% of the maximum light intensity. 

10 

1 1 Just before the storm, the MHD E5D simulation predicts at 00:00 UT a plasmasphere 

12 extending up to 6 Re. Since the dawn-dusk component of the two other convection 

13 electric fields (VMSC and Weimer) are stronger than that of the E5D model, the LCS is 

14 located closer to the Earth for these two other models: i.e. around 4 Re for the MHD 

15 VSMC, and even closer for the MHD W simulation. When the interchange mechanism is 

16 assumed to peel off the plasmasphere, the plasmapause at midnight LT is located around 

17 4.5 Re with the E5D model; this is closer to Earth than the LCS simulation for this same 

18 E-field model. However, the plasmapause predicted by the IC E5D simulation is located 

19 at larger equatorial distances than the LCS for the two other E-field models. But for a 

20 same E-field model, the minimum equatorial distance of the ZPF surface is generally 

21 smaller than that of the LCS. 

22 

23 Note that the plasmasphere is almost circular at that time which corresponds to the initial 

24 time to in all different simulations. Although the level of geomagnetic activity has been 

25 nearly constant and quite low during the previous 24 h, a bulge is present at t 0 in the dusk 

26 sector in the MHD VSMC simulation. However, according to the EUV observations the 

27 plasmasphere is usually circular after prolonged quiet periods. Plasma tails or plumes are 

28 formed only during disturbed periods [Pierrard and Cabrera, 2005]. Unfortunately, no 

29 observations of EUV/IMAGE are available at 00:00 UT 28 October 2001, since the 

30 satellite was not close to its apogee at that initial time to. 

31 
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1 At 04:30 UT as illustrated in Fig. 10, the geomagnetic activity level Kp increases 

2 significantly and reaches a maximum Kp = T. This increase of Kp is associated to a 

3 southward turning of the interplanetary magnetic field Bz, as well as a decrease of Dst 

4 index. The latter reaches a minimum of- 157 nT at 1 1 :00 UT. By that time, the 

5 convection electric field intensity has increased in all models so that the corresponding 

6 LCS-24h are less extended. Only the innermost flux tubes have a circular trajectory, 

7 while the more distant closed streamlines are highly skewed with a maximum radial 

8 distance in the post-noon sector. Beyond the LCS, the plasma elements are lost to the 

9 magnetosheath and are not any longer tracked in our MHD simulations. According to this 

10 MHD theory for the formation of the plasmapause, a more or less sharp knee is formed 

1 1 along the LCS-24h, wherein the streamlines have remained closed for at least 24 hours. 

12 The LCS shrinks to about 3 Re in the night sector for all MHD simulations. A bulge is 

13 formed in the afternoon sector for all the simulations as a result of the enhancement of the 

14 dawn-dusk E-field component. Indeed, this enhancement of E-field in the dusk sector 

1 5 produces a sunward surge of plasma in this MLT sector. Note, however, the quite 

16 different shapes of the equatorial cross section of the LCS predicted by the various E- 

17 field models. The bulge in the LCS is the precursor of the plume that usually develops 

18 during geomagnetic storms and substorms. 

19 

20 The interchange mechanism shown in the fourth panels produces much more irregular 

21 plasmapause shapes than the MHD simulations. From the fourth panel in Fig. 10, it can 

22 be seen that the ZPF surface where the plasmasphere is peeled off in the postmidnight 

23 sector corresponds to an equatorial distance slightly beyond 2 Re. Despite the rather 

24 limited quality of the observations of IMAGE at the time of 04:30 UT, it can be seen that 

25 the plasmapause is located close to the Earth during this storm. 

26 

27 At 08:00 UT, the Kp index is still high: Kp = 6'. As illustrated in Fig. 1 1 , a plume is now 

28 clearly developed in all the simulations, with and without interchange. The plumes are 

29 located in the afternoon sector but have rather different shapes and experience quite 

30 different development depending on the adopted E-field model. The observation of 
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1 EUV/IMAGE at 08:03 UT is very contaminated, but a plume can clearly be identified in 

2 the afternoon sector. 

3 With MHD VSMC and Weimer simulations, the plasmasphere is still smaller than with 

4 MHD E5D, due to their stronger dawn-dusk E-field intensity which implies a stagnation 

5 point rather close to Earth. The simulation of interchange with the E5D model, gives a 

6 position of the plasmapause which is rather close to that of the LCS of MHD E5D 

7 simulation. 

8 

9 By 12:00 UT, the plume observed in the simulations has slightly rotated eastward as can 

10 be seen in Fig. 12. No EUV observations are unfortunately available at that time. At 

1 1 18:05 and 19:58 UT, EUV/IMAGE provides nice observations, and a plume is clearly 

12 visible in the aftemoon/dusk region as illustrated in Figs. 13 and 14. The results of the 

13 simulations are also shown for 18:00 and 20:00 UT. The longitudinal extent and MLT 

14 position of the plume is again slightly different in the different models. Note that by 

1 5 reducing arbitrarily the co-rotation velocity, a better agreement would be obtained. This 

1 6 kind of forcing has been avoided in the present work in order not to distort the original E- 

1 7 field models which have been derived without reducing the co-rotation electric field 

1 8 intensity. Note nevertheless that the velocity of the plasmasphere is observed to be 

19 generally slower than co-rotation [Burch et al., 2004], what explains the faster rotation of 

20 the plume in the different simulations compared to the EUV observations. 

21 

22 Similar plumes are often observed by EUV/IMAGE and CLUSTER after a significant 

23 increase of the geomagnetic activity level [Darrouzet et al., 2006a and b]. They are 

24 formed in the afternoon sector and then rotate eastward with the core of the 

25 plasmasphere. Other examples of IMAGE observations during quiet and disturbed 

26 periods that have been compared to results of simulations based on the interchange 

27 mechanism (IC E5D) have been presented by Pierrard and Cabrera [2005, 2006] and 

28 Pierrard [2006]. 

29 

30 4.2 Magnetic storm of 17 April 2002 


20 



1 Fig. 15 shows the results of similar simulations obtained during the geomagnetic storm of 

2 17-4-2002. This case is particularly interesting since this is the single event for which the 

3 effect of several electric field models on the position of the plasmapause has been also 

4 determined by Liemohn et al. [2004], The E-field models tested in this paper were (i) 

5 their modified E5D model rescaled to have a larger intensity than Mcllwain’s original 

6 version, (ii) the Weimer model, and (iii) a self consistent E-field determined by the 

7 authors. Thus it is interesting to compare the results of the present MHD simulations 

8 based on VSMC as well as on the original E5D model version with the results found by 

9 Liemohn et al. [2004] for the same event on 1 7-4-2002 with other E-fields models. 

10 Moreover, we also show the effect of the interchange mechanism on the positions, shape 

1 1 and evolution of the plasmapause during this magnetic storm. 

12 

13 Fig. 15 shows the results at 21:00 UT on 17-4-2002 obtained with the MHD simulations 

14 respectively for E5D, VSMC and Weimer models (three upper panels), as well as with 

15 interchange mechanism based on the E5D model (left bottom panel). A plume is again 

16 formed with all the electric field models. Its development is associated to the increase of 

17 the geomagnetic activity level up to Kp — T . This increase of Kp is also associated with a 

1 8 southward turning of the interplanetary magnetic field, like during all geomagnetic 

19 storms. During the main phase, Dsit decreases reaching a minimum of -105 nT on 17-4- 

20 2002 and even a lower value the day after. The plume is located approximately in the 

21 same MLT sector as in the observations of EUV. The MHD E5D simulation shows a 

22 plasmapause quite similar to that obtained with the rescaled E5D model used by Liemohn 

23 et al. The position of the plasmapause found in Liemohn et al. MHD simulation was in 

24 good agreement with the observations of EUV in the nightside, but it was too far from the 

25 Earth on the dayside. We obtain the same characteristics with the MHD E5D simulation. 

26 The MHD VSMC simulation corresponds better to the plasmapause observations on the 

27 dayside, but on the nightside, it is too close from the Earth. The MHD W simulation 

28 based on Weimer model gives also a plasmapause closer to the Earth than what is 

29 observed, especially on the nightside; these results are in agreement with the simulation 

30 published by Liemohn et al. where the Weimer convection E-field is also used. 

31 
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1 The simulation of the interchange mechanism with the E5D model gives results closer to 

2 the EUV observations than the MHD simulations with this same E5D electric field 

3 model. This does not imply, of course, that there is no need for better E-field model than 

4 E5D, nor is it a pleading for the mechanism of interchange, since any physical theory is 

5 necessarily an approximation and is therefore perfectible. 

6 

7 While in general, interchange simulations with E5D model predict results that are in good 

8 agreement with the EUV observations during substorm events [Pierrard and Cabrera, 

9 2005], we have seen here that such an agreement is more precarious during the 

10 geomagnetic storms examined in this paper. We attribute the lack of satisfactory results 

1 1 during geomagnetic storms to the lack of characterization of the E5D electric field model 

12 under this sort of geomagnetic disturbances. Indeed, this empirical E-field model was not 

13 designed to reproduce the magnetospheric convection electric field during geomagnetic 

14 storms, but only to model the dispersion of electrons and protons accelerated or injected 

15 into the magnetosphere following substorm events [Mcllwain, 1986]. This lays the need 

16 to update and improve currently available empirical electric field models not only at 

17 ionospheric level but also at high altitudes in the Earth magnetosphere. Furthermore, 

1 8 nobody can argue that future alternative theories for the formation of the plasmapause 

19 will not be able to simulate the dynamics of the plasmasphere more accurately and offer 

20 predictions in closer agreement with observational reality. 

21 

22 5. Conclusions 

23 

24 After brief descriptions of different electric fields models (: the Volland-Stem-Maynard- 

25 Chen model VSMC, Mcllwain’s E5D model and the Weimer model), we first determined 

26 how the distribution of the equatorial equipotentials corresponding to these 

27 magnetospheric convection E-field models change during the geomagnetic storm of 28 

28 October 200 1 . Their equipotentials in co-rotating and non-corotating frames of reference 

29 depend on the level geomagnetic activity; they are thus changing as a function of 

30 universal time since the geomagnetic activity indexes Kp and Dst are varying hour after 

31 hour. 
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1 The radial and local time distributions of these equipotential contours have been 

2 displayed in Figs. 1, 2, 4, 5, 6, and 7 for a set of times for which EUV observations are 

3 sometimes available. The E-field models are used in association with magnetic field 

4 models (VSMC/dipole, E5D/M2 model, and’Weimer/Tsyganenko respectively), to 

5 calculate the convection velocity of cold plasma or the drift velocity of zero energy 

6 charged particles. The dipolar and M2 magnetic field models had been employed to 

7 design, respectively, the VSMC and E5D models from in-situ satellite observations. In 

8 addition to their portability, these simple quasi-stationary E-field models have an 

9 expedient advantage: they depend only on one single parameter, the three-hourly Kp 

10 geomagnetic activity index. This makes them user friendly and easy to implement in any 

1 1 numerical codes. Weimer E-field is more sophisticated: it depends on several solar wind 

12 parameters and varies with a smaller time resolution. Nevertheless, it is based on low 

13 altitudes and high latitudes observations collected far from the equatorial region where 

14 the plasmapause is formed. 

15 Although all these empirical electric field models are based on observations (respectively, 

16 on OGO-3 & 5 satellites for VSMG, AST-5 & 6 satellites for E5D and low altitude 

17 ionospheric convection velocity measurements at high latitudes for Weimer), they occur 

18 to be quite different from each other. Note that these models are averages and 

19 approximations of the actual field distributions at any particular instant of time, and at 

20 any particular place in the magnetosphere. Some of them were not designed to model 

21 magnetospheric electric fields in cases of rapidly changing magnetic field intensities: e.g. 

22 during geomagnetic storms when induction electric fields are generated on top of the 

23 electrostatic component approximated by these E-field models. For instance, the E5D 

24 model was essentially built to represent the E-field distribution immediately after a 

25 substorm injection event; furthermore, it was developed to represent the E-field in the 

26 region of geosynchronous orbit, only when the level of geomagnetic activity, denoted by 

27 the value of Kp, remains nearly constant and smaller than 6. Mcllwain’s electric field 

28 E5D was not developed under such circumstances as geomagnetic storms. But, taking 

29 into account these warnings and restrictions, it is interesting to note that the plasmapause 

30 position obtained with this E-field model and the IC mechanism is in better agreement 
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1 with the EUV observations during the geomagnetic storm of 1 7 April 2002 than the 

2 positions obtained with the other E-fields models and MHD simulations. 

3 

4 Two mechanisms of plasmapause formation wefe also compared and confirm that the 

5 last-closed-streamline (LCS) of the MHD simulation is located beyond the plasmapause 

6 predicted by the interchange mechanism for the geomagnetic storm of October 28, 2001 : 

7 the latter being thus generally closer to the Earth. A similar comparison with results 

8 obtained for the geomagnetic storm of April 17, 2002 leads to the same conclusions. 

9 

10 The choice of the E-field model is crucial in the results of the simulations. Plumes 

1 1 develop during geomagnetic storms and substorms with all the various E-field models, 

12 but quite different shapes of the equatorial cross section of the LCS are predicted. It is 

13 quite clear that more detailed magnetospheric E-field models with higher time resolutions 

14 should be developed. In such future E-field models the distribution of the equipotential 

15 surfaces should possibly be desynchronized: their evolution in the night side should not 

16 necessarily be synchronized with that in the dayside or at any other MLTs. Finally, it 

17 would be useful to model also the inductive electric field component in order to model 

1 8 more properly what happens during geomagnetic storms. 

19 Despite the better scores of the E5D electric field model and of the interchange scenario 

20 for the formation of the plasmapause found in many case studies, it must be admitted, 

21 however, that the equatorial cross section of the plasmapause as determined from the 

22 EUV observations during both geomagnetic storms, are not very well reproduced by none 

23 of the simulations presented above, not even when the Weimer electric field model is 

24 used instead of the E5D model. This leads us to conclude that none of the electrostatic 

25 field models is fully adapted to model the actual magnetospheric E-field distribution 

26 during geomagnetic storms (at least for those selected in the present study). 

27 

28 Since previous simulations by Pierrard and Lemaire [2004] as well as by Pierrard and 

29 Cabrera [2005] based on the interchange mechanism and the E5D electric field model 

30 have shown that overall shapes of the plasmapause and its evolution does rather well 

31 explain the formation of plasmatails or plumes as well as shoulders following to 
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1 substorm events, it may be speculated that during geomagnetic storms characterized by 

2 larger Dst variations, the time dependent electric field distribution has an induced 

3 component associated to the increase of southward magnetic field which is generated by 

4 the King Current during the main phase. A toroidal induced electric field of smaller 

5 intensity and of opposite direction is also expected during the recovery phase of 

6 geomagnetic storms. These toroidal induced electric fields cannot be represented as the 

7 gradient of a scalar electrostatic potential as most empirical magnetospheric E-field 

8 models used above. The effect of such additional time dependent and non curlfree electric 

9 field distribution on the plasmasphere and on the formation of the plasmapause, has not 

1 0 yet been evaluated in detail. 

11 It is suggested here that such induced electric fields generated during geomagnetic storms 

12 are responsible for the lack of satisfactory agreement between the simulations presented 

13 in this study based on curlfree electric field models. We suggest that this effect should be 

14 examined in the future and included in forthcoming and more comprehensive theories 

1 5 for the formation of the plasmapause not only during substorm events but also even 

1 6 during geomagnetic storms with large Dst variations. 

17 

1 8 Therefore, this study points out the need to develop higher time resolution empirical 

19 models for the magnetospheric electrostatic field distribution like those developed for the 

20 geomagnetic field. It urges to take into account the effect of induced electric fields 

21 generated in the magnetosphere during geomagnetic storms with large Dst variations. It is 

22 only when such more detailed and comprehensive E-field models will be available in 

23 association with the time dependent empirical B-field models, that one might expect to 

24 test quantitatively and more definitely any of the existing and future theories for the 

25 formation of the plasmapause by comparing their theoretical predictions to observations 

26 like those of EUV. 

27 
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32 Figure captions 

33 

34 Figure 1 : Equipotential contours of the Kp-dependent convection electric field Volland- 

35 Stem on 28 October 2001 (without co-rotation electric field). 

36 

37 Fig. 2: Equipotential map of the Kp-dependent convection electric field E5D on 28 

38 October 2001 (without co-rotation electric field). 

39 

40 Fig. 3: Magnetic field model M2. 

41 

42 Fig. 4: Equipotential map of the Weimer convection electric field dependent on solar 

43 wind conditions on 28 October 2001 (without co-rotation electric field). 

44 

45 Fig. 5: Equipotential map of the total electric field Volland-Stem on 28 October 2001 . 

46 Same as Fig. 1 but with co-rotation. 

47 

48 Fig. 6: Equipotential map of the total electric field E5D on 28 October 2001. Same as 

49 Fig. 2 but with co-rotation. 

50 

5 1 Fig. 7 : Equipotential map of the Weimer total electric field dependent on 28 October 

52 2001 . Same as Fig. 4 but with co-rotation. 
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1 

2 Fig. 8: Variation of the solar wind density of protons, the solar wind velocity and the Bx 

3 and By component of the magnetic field. 

4 

O ji t;;.J 

6 Fig. 9: Initial plasmasphere determined with the MHD simulations using the E5D E-field 

7 model (upper left panel), MFID with Volland-Stem model (upper middle panel), MHD 

8 with Weimer model (upper right panel) and the plasmapause position predicted by 

9 instability mechanism with E5D (bottom left panel) on 28 October 2001 at 0:00 UT. All 

10 these simulations started at 0:00 UT, 27 October 2001. The upper panels show the values 

11 of the Northward component of the interplanetary magnetic field (Bz), the value of the 

12 ring current magnetic field (Dst in nT), and the values of the geomagnetic index Kp 

13 during 3 days between 27 October 2001 and the end of 29 October 2001 . No observations 

1 4 from EUV/IMAGE are available at that time. 

15 

16 Fig. 10: Results of the MHD E5D simulations (upper left panel), MHD VSMC (upper 

17 middle panel), MHD W (upper right panel), IG E5D (bottom left panel) at 4:30 UT, 28 

1 8 October 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 

19 4:39 UT is illustrated in the bottom right panel. 

20 

21 Fig. 11: Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 

22 middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 8:00 UT, 28 

23 October 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 

24 8:03 UT is illustrated in the bottom right panel. 

25 

26 Fig. 12: Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 

27 middle panel), MHD W (upper right panel), IC E5D (bottom panel) at 12:00 UT, 28 

28 October 2001. No EUV/IMAGE observation is available at that time. 

29 

30 Fig. 13: Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 

3 1 middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 1 8 :00 UT, 28 

32 October 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 

33 1 8:05 UT is illustrated in the bottom right panel. 

34 

35 Fig. 14: Results of MHD E5D simulations (upper left panel), MHD VSMC (upper middle 

36 panel), MHD W (upper right panel), IC E5D (bottom left panel) at 20:00 UT, 28 October 

37 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 1 9:58 UT 

38 is illustrated in the bottom right panel. 

39 

40 Figure 15: Results of MHD E5D simulations (upper left panel), MHD VSMC (upper 

41 middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 21 :00 UT, 

42 after the magnetic storm of 17 April 2002. The EUV/IMAGE observation of the 

43 equatorial plasmapause position at 21 :07 UT is illustrated in the bottom right panel. 

44 
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Figure 1 : Equipotential contours of the Kp-dependent convection electric field Volland- 
Stem on 28 October 2001 (without co-rotation electric field). 
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Fig. 2: Equipotential map of the Kp-dependent convection electric field E5D on 28 
October 2001 (without co-rotation electric field). 
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1 Fig. 4: Equipotential map of the Weimer convection electric field dependent on solar 

2 wind conditions on 28 October 2001 (without co-rotation electric field). 
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3 Fig. 5: Equipotential map of the total electric field Volland-Stem on 28 October 2001 

4 Same as Fig. 1 but with co-rotation. 
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Fig. 6: Equipotential map of the total electric field E5D on 28 October 2001. Same as 
Fig. 2 but with co-rotation. 
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1 Fig. 7: Equipotential map of the Weimer total electric field dependent on 28 October 

2 2001 . Same as Fig. 4 but with co-rotation. 
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1 

2 Fig. 9: Initial plasmasphere determined with the MHD simulations using the E5D E-field 

3 model (upper left panel), MHD with Volland-Stem model (upper middle panel), MHD 

4 with Weimer model (upper right panel) and the plasmapause position predicted by 

5 instability mechanism with E5D (bottom left panel) on 28 October 2001 at 0:00 UT. All 

6 these simulations started at 0:00 UT, 27 October 2001 . The upper panels show the values 

7 of the Northward component of the interplanetary magnetic field (Bz), the value of the 

8 ring current magnetic field (Dst in nT), and the values of the geomagnetic index Kp 

9 during 3 days between 27 October 200 1 and the end of 29 October 2001 . No observations 

1 0 from EUV/IMAGE are available at that time. 
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Fig. 10: Results of the MHD E5D simulations (upper left panel), MHD VSMC (upper 
middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 4:30 UT, 28 
October 2001. The EUV/IMAGE observation of the equatorial plasmapause position at 
4:39 UT is illustrated in the bottom right panel. 




39 


1 

2 Fig. 1 1 : Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 

3 middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 8:00 UT, 28 

4 October 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 

5;, 8:03 UT is illustrated in the bottom right panel. 
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Fig. 12: Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 
middle panel), MHD W (upper right panel), IC E5D (bottom panel) at 12:00 UT, 28 
October 2001 . No EUV/IMAGE observation is available at that time. & 
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1 

2 Fig. 13: Results of the MHD E5D simulation (upper left panel), MHD VSMC (upper 

3 middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 1 8:00 UT, 28 

4 October 2001 . The EUV/IMAGE observation of the equatorial plasmapause position at 

5 18:05 UT is illustrated in the bottom right panel. ,, 

6 
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1 Fig. 14: Results of MHD E5D simulations (upper left panel), MHD VSMC (upper middle 

2 panel), MHD W (upper right panel), IC E5D (bottom left panel) at 20:00 UT, 28 October 

3 2001. The EUV/IMAGE observation of the equatorial plasmapause position at 19:58 UT 

4 is illustrated in the bottom right panel. 
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Figure 15: Results of MHD E5D simulations (upper left panel), MHD VSMC (upper 
middle panel), MHD W (upper right panel), IC E5D (bottom left panel) at 21 :00 UT, 
after the magnetic storm of 17 April 2002. The EUV/IMAGE observation of the 
equatorial plasmapause position at 21 :07 UT is illustrated in the bottom right panel. 
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